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Abstract 

Using an analogy with the well-known double-slit experiment, we show that the 
standard phase of neutrino oscillations is correct, refuting recent claims of a factor of 
two correction. We also improve the wave packet treatment of neutrino oscillations 
taking into account explicitly the finite coherence time of the detection process. 



The experimental and theoretical investigation of neutrino oscillations is presently one of 
the most active fields of research in high-energy physics. For the interpretation of experimen- 
tal data it is important to have a clear and correct theoretical formulation of the neutrino 
oscillation mechanism^} In particular, the phase &kj of neutrino oscillations in vacuum due 
to the interference of the contributions of two massive neutrinos Uk and Uj with masses rrik 
and rrij, respectively, depends on the mass-squared differences AmjL = ml — mj. A correct 
theoretical expression of the oscillation phase &kj is important in order to extract from the 
experimental data correct information on the value of AmjL. 

The standard expression for the oscillation phase $>kj i n the relativistic approximation is 

where E is the neutrino energy and L is the distance between the neutrino source and 
detector. 

It has been argued by some authors [|3],[|] that the expression ([I]) is wrong by a factor of two 
on the basis of the following reasoning. Neutrino experiments measure neutrino oscillations 
as a function of the source-detector distance L. Different massive neutrinos propagate with 
different velocities 

Pk / v 
V k = TT> (2) 

&k 

where E k and p k are, respectively, the energy and momentum of the neutrino with mass m&, 
related by the relativistic dispersion relation 

El=pl + m 2 k . (3) 

Hence, the phases of the different massive neutrinos wave functions after a propagation 
distance L should take into account the different times of propagation of different massive 
neutrinos: 

$ k =p k L- E k t k . (4) 
The different propagation times are given by 

t k = ± = *L, (5) 
v k Pk 

which lead, in the relativistic approximation, to the phase difference 

~ Am?, L 
A$ kj = $ fc - = j*- . (6) 

This phase difference is twice of the standard one in Eq. ([[]). A similar disagreement by a 
factor of two has been claimed to exist in the case of kaon oscillations in Refs. 

Let us notice that in Eq. (^) we have considered the possibility of different energies and 
momenta for different massive neutrino wave functions. Indeed, it has been shown in Ref. || 



1 See Ref. § for a review of the theory of neutrino oscillations and Ref. Q for an exhaustive list of 
references. 
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that Lorentz invariance implies that in general in oscillation experiments different massive 
neutrinos have different energies and different momenta. The energy E in Eqs. (H) and (§) 
is the energy of the massive neutrinos in the massless limit, i.e. neglecting the differences of 
energy due to the masses. 

The authors of Refs. [p|-|TT|| claimed that a correct way to obtain the standard oscillation 



phase is to assume the same energy for the different massive neutrino wave functions. This is 
an unphysical assumption, as discussed in Ref. M. Moreover, as already noticed in Ref. [Q, 
it is not true that the disagreement of a factor of two disappears assuming the same energy 
for the different massive neutrino wave functions, as clearly shown by the above derivation 
of Eq. (^|), in which the energies of the different massive neutrino wave functions could have 
been taken to be equal. Indeed, even if the different massive neutrino wave functions have 
the same energy, the time contribution — Et k + Etj to the phase difference AQ k j does not 
disappear, because t k 7^ tj. This contribution has been missed in Refs. [PHTTH. 



In order to test the validity of the reasoning of the authors of Refs. |3|,f|], let us apply 
their method to the well-known double-slit interference experiment depicted in Fig. [I]. The 
particles emitted by the source S and detected could be photons or electrons, or others. 
The screen A has two holes through which the particles can reach the screen B on which a 
detector registers the arrival of the particles at a distance x from the center of the screen. 

In the standard approach, the phases of the two waves at x at the time t are given by 

® k =pr k -Et = p(r A + r Bk )-Et, (7) 
for k = 1,2. Hence, the phase difference is 

A$E$ 2 -$ 1= pAr, (8) 

with 

2 a x 

Ar = r 2 — r x = r B 2 — r B i — — r~ j (9) 

a 

where the approximation holds for x,a -C d. Since p = 2tt/X, where A is the wavelength, 
the maxima of interference are given by the usual well known formula x = n\d/2a, with 
n = 0, 1, . . ., which has been confirmed by many experiments without any doubt. 

Let us apply now the reasoning of the authors of Refs. [|, |]. If one takes into account 
the velocity of the particles, the two paths have different propagation times 

** = -, (10) 

v 



where v is the velocity, 



In this approach, the phases of the two waves at x are given by 

$ k =pr k -Et k . (12) 
The phase difference at x turns out to be 

~ ~ ~ / E\ . am 2 , . 

A$ = $ 2 -$i= \P Ar~-2-x — , (13) 

\ v I dp 



2 



where the approximation holds for x,a <C d. This phase difference is very different from the 
correct one in Eq. (J|), and it even vanishes for massless particles (as photons). 
Therefore, the reasoning presented in Refs. is wrong. 

The mistake in the claim formulated in Refs. is due to a wrong use of the group 
velocity in the phase, which depends on the phase velocity. The group velocity (given in 
Eq. (0) for neutrinos and in Eq. (|TT| ) for the double-slit experiment) is the velocity of the 
factor that modulates the amplitude of a wave packet describing a localized particle. In the 
double-slit experiment there are two wave packets which propagate along the two paths in 
Fig. [I]. The envelopes of these wave packets take different times to cover the two different 
distances from the source S to the point x on the screen. But this has no effect on the phases. 
Only the amplitude of the final wave function is determined by the modulating factors of 
the amplitudes of the two wave packets. The different arrival times of the envelopes of the 
wave packets reduces the overlap of the two wave packets, leading to a decoherence effect. 

Let us illustrate these concepts in a one-dimensional formalism using a Gaussian momen- 
tum distribution with width o p , centered at the momentum p: 



2n o, 



-1/2 



exp 



(p' - pY 

4o 2 
p 



(14) 



For Op « p, the corresponding wave packet in space-time is 



ip{r,t) 



1/2 

2n o r ) exp 



ipr — iEt — 



(r — vty 

4o-2 



(15) 



where o r = l/2o p , E = yp 2 + m 2 , and v = p/E is the velocity of the envelope 



-1/2 



exp 



(r — vty 

4o 2 



(16) 



of the wave packet. The real part of the wave packet ([15|) at t — is depicted by the solid 
line in Fig. |2|, with the envelope flT6|) represented by the dashed line. It is important to 
understand that the wave packet in the approximation ([15]) is a monochromatic wave whose 
amplitude is modulated by the envelope factor (|T6|) , that has no effect on the phase factor 



exp [ipr — iEt] 



(17) 



of the wave packet. 

The amplitudes of each of the two wave packets traveling the two paths in the double-slit 
experiment in Fig. [l| at the distance x from the center of screen B at the time t is given by 



ip{r k ,t) = -j= (y2na r 



-1/2 



exp 



iprk — i Et — 



(r k - vty 

4o 2 



for k = 1,2, with r k = ta + tbu- The probability to detect the particle at the point x at the 
time t is 



P(x,t) oc \*l>{r x ,t)+<il>{r 2 ,t)f 
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oc 



2 V27T0V 



cxp 



(ri — viy 



2a 2 r 



exp 



(r 2 - vty 



+2 cos(pAr) exp 



2a 2 r 

\2 , / ,\2 



(r x - ft) + (r 2 - vty 
4^! 



(19) 



where we have omitted an appropriate normalization factor, that must be inserted in order 
to normalize the probability over the screen. 

From Eq. ([19]) it is clear that the probability to detect the particle at x is appreciable only 
when one of the two wave packets cover the detector. Interference of the two wave functions 
is possible if both wave packets cover the detector simultaneously. However, usually the 
detector operates for a long time without recording the precise instant of detection. In this 
case, the probability of detection at x is given by the time average of P(x,t) in Eq. ([TP]): 



P(x) oc — >{ 1 + cos(pAr) exp 



Ar 5 



(20) 



The interference term gives an appreciable contribution only if the separation Ar of the two 
wave packets is smaller than their width, i. e. if there is an overlap of the two wave packets 
at the detector. It is clear that the different propagation times of the envelopes of the two 
wave packets has no effect on the phase of the interference term, that is equal to the phase 
in Eq. (g), obtained in the plane wave approximation. This is what we wanted to show with 
this simple example. 

An interesting complication to the problem stems from the fact, noted in Ref. [12], that 
different wave packets can interfere even if they do not overlap at the detector if the coherence 
time of the detection process is larger than the time separation between the arrival of the 
wave packets. In order to take into account this effect in the double-slit experiment, we 
have to introduce an appropriate wave function for the detection process. Assuming that 
the detection process has a coherence time a t r> and the detection occurs at the time fro, we 
describe the detection process with the Gaussian wave packet 



^d(Md) 



27T cr tD 



-1/2 



cxp 



-iE(t-t D ) 



(t - t D y 

4^ 2 D 



(21) 



The frequency of the detection process is the same as that of the absorbed wave function, 
because the incoming wave function excites the corresponding degree of freedom of the 
detection process. The amplitude of detection of each of the two incoming wave packets is 
given by the overlap of the incoming wave packet (pj) and the detection process wave packet 



Performing the integration over t, we obtain 



A(r k ,t D ) oc — f= (V2~ko 



A(r k ,t B ) oc / dtip^(t,t D ) ip(r k ,t) . 



ipr k -iEt D - 



-1/2 



exp 



(r k - vt D ) 

4(T 2 



(22) 



(23) 
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with 



a- = a 2 + v 2 a 2 D . (24) 

One can see that the time allowed for a coherent absorption of the incoming wave packet is 
increased by the coherence time of the detection process, that dominates if v a t D ^> ay. 
The probability to detect the particle at x at the time is now given by 



P(x,t D ) oc \A(r u t D ) + A(r 2 ,t D )f 



1 



2V2na 



exp 



(ri - vt D y 

2a 2 



exp 



(r 2 -i;t D )' 
2a 2 



+2 cos(pAr) exp 



(ri - ^tp) 2 + (r 2 - ftp)' 

4(T 2 



(25) 



This equation has the same structure as Eq. (0), with the width oy of the wave packets 
replaced by a. Obviously, the average over the unmeasured detection time tp leads to an 
expression for the probability to detect the particle at x given by Eq. ([HJ) with oy replaced 
by a, which takes into account also the coherence time of the detection process. 

Therefore, we have seen that the coherence time of the detection process, which allows a 
coherent absorption of wave packets arriving at the detector at different times, do not have 
any effect on the phase of the interference between different wave functions, refuting the 
arguments presented in Refs. [|3],[§]. The resulting prescription in calculations performed in 
the plane wave approximation (i.e neglecting the wave packet character of wave functions 
that describe localized particles) is to calculate the interference of different wave functions 



at the same time and the same space point [13,14]. 

The physical explanation of the independence of the phase of the interference between 
different wave functions from their arrival time at the detector has been presented in Ref. ||15|| . 
It consists in taking into account that there is interference only if the different wave functions 
are detected coherently. Coherence means that there is a precise phase difference of the 
detection process between the arrival times of the different wave functions that must be 
taken into account. When the second wave packet arrives, the detection process is already 
excited with a frequency due to its interaction with the first wave packet and its phase is 
determined by the time difference between the arrivals of the two wave packets. It is easy 



to see |Lo] that the phase of the detection process cancels exactly the phase difference of the 
two wave functions due to the different detecting times, leading to the practical prescription 
to calculate the interference of different wave functions at the same time and the same space 
point. 

Let us now apply the wave packet formalism to neutrino oscillations, following the lines 
presented in Refs. [|1(| [17]]. The Gaussian wave packets that describe massive neutrinos in 
one spatial dimension x along the source-detector direction are 



ip k (x,t) 



27TCT 



xP 



'1/2 



exp 



ip k x -iE k t — 



(x - v k ty 

Kp 



(26) 



where a x p is the width of the wave packets, that is determined by the production process, E k 
and pk are, respectively, the energy and momentum of the neutrino with mass m^, related 
by the relativistic dispersion relation (H), and is the group velocity given in Eq. (0). 
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Taking into account the possibility that the detection process has a finite coherence time 
interval <r t D and a finite spatial coherence width a xD , we describe the detection process of 
the massive neutrino v k with the Gaussian wave packet 



ipT) k (x,t,L,T) oc exp 



ip k {x - L) -i E k (t - T) 



{x-Lf (t-Tf 



Kb 



Kb 



(27) 



where L is the source-detector distance and T is the time elapsed between neutrino produc- 
tion and detection. Notice that with respect to the wave packet (PIP, in Eq. ( p7| ) we have 
considered also a finite spatial coherence width of the detection process that was neglected 
in the discussion of the double-slit experiment. On the other hand, with respect to the wave 
packet treatment of neutrino oscillations presented in Ref. [17 we have added the explicit 
terms that take into account the finite coherence time interval a t r> of the detection process. 

The frequency and wave number of the detection process are the same as those of the 
detected massive neutrino wave function, which excites the corresponding degree of freedom 
of the detection process. In the framework of quantum field theory the detection process is 
described in terms of the leptonic weak charged current 



a=e,iM,r k 



^a(x) 7 M (I-7 5 ) U ak u k (x) 



(2? 



where U is the lepton mixing matrix (see [jTB|-|2T||). Only the field v k {x) is excited by the 
arrival of the wave function of the corresponding massive neutrino u k , and the frequency and 
wave number of the field excitations are the same as those of the incoming wave function. 

The amplitude of detection of each massive neutrino wave packet is given by the overlap 
of the incoming neutrino wave packet with the detection process wave packet: 



A k (L,T) oc J dxdt 



(29) 



Performing the integration over x and t we obtain 



A k (L, T) oc exp 



ip k L -%E k T - 



[L - v k Tf 
Kk 



with 



°4 = °lv + °Ib + V l a tD 



(30) 



(31) 



This result is simple and important. It shows clearly that the coherent absorption of each 
massive neutrino wave function occurs in a space-time interval of width a xk around the 
coordinates L, T of the detection process. The size of this space-time interval is determined 
by the width of the neutrino wave packet and the coherence spatial and temporal widths of 
the detection process. The largest width dominates. The contribution of the temporal width 
of the detection process is weighted by the velocity of the incoming neutrino wave packet 
for obvious reasons. Since this velocity depends on the mass of the neutrinos, the width a xk 
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depends on the index k labeling massive neutrinos. However, for relativistic neutrinos the 
contribution of neutrino mass to a xk is negligible and we can safely approximate 



2 2 



, _2 , 2 



(32) 



Notice that, on the other hand, the contribution of neutrino mass to the other terms in 
Eq. ( j30|) cannot be neglected because it is amplified by the macroscopic quantities L and T. 
In the following we use the relativistic approximation, which implies that 



Eh 



E + ^ 
^ 2E 



Vk 



E 



W 2E 



v k 



m k 
2E 2 



(33) 



where E is the neutrino energy in the limit of zero mass. The first order correction to 
the momentum and energy due to neutrino mass m k is proportional to m\ because of the 
relativistic dispersion relation ([|), and must be divided by the energy E for dimensional 
reasons. The coefficient £ depends on the production process, but we will see that it has no 
effect on the phase of neutrino oscillations, as already shown in Refs. [I5-I7|. 



The probability of flavor transitions is given by (see [18-21]) 



P Va ^ u JL,T) oc 



J2KkML,T)U, 



i3k 



(34) 



where U is the lepton mixing matrix. Performing the average of P Ua ^ U/3 (L,T) over the 
unmeasured neutrino propagation time T, we finally obtain 



P Va ^ V0 (L) = J2\Uak | 2 1 U pk 1 2 + 2 Re Kk Up k U aj U* 0j F k j exp 

k k>j 



where 



are the coherence lengths and 



■ coh 



AV2E 2 



Am 2 kj L ( L 



2E 



T coh 
L kj 



(35) 
(36) 



F kj = exp 



-2^i 



2 t-2 / ° & 



T OSC 
L kj 



with the oscillation lengths 



T OSC 



AttE 
\ Am lj\ 



(37) 



(38) 



The coherence length L£° h is the distance beyond which the contributions of the massive 
neutrinos v k and Uj do not interfere any more because the separation between their wave 
packets is larger than the size of the wave packets and the spatial and temporal coherence 
widths of the detection process. The coefficient F k j suppresses the interference of v k and Vj 
if a x > L^ c , i.e. if the production or the detection process is not localized in a space-time 
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region much smaller than the oscillation length. This obvious constraint for the observation 



of neutrino oscillations was discussed for the first time in Ref. |22|] and is satisfied in all 
experiments. Therefore, in practice one can safely approximate F^j — 1. 

From Eq. (j35|) one can see that a proper wave packet treatment of neutrino oscillations 
that takes into account the different propagation times of massive neutrino wave packets leads 
to the standard oscillation phase (|l|), refuting the disagreement of a factor of two claimed in 
Refs. ||,£| . Neutrino oscillations are due to the different phase velocities of different massive 
neutrinos, which produces interference. This interference is obviously the same whether it 
is calculated in the plane wave approach or in the wave packet treatment. The different 
conclusion reached in Refs. |3|,f3| is due to a wrong use of the group velocity in the phase. 

In conclusion, using an analogy with the well-known double-slit experiment we have 
shown that the standard phase ([!]) of neutrino oscillations is correct, refuting the claim of 
a factor of two correction presented in Refs. P,ffl. We have also improved the wave packet 



treatment of neutrino oscillations presented in Refs. |T6|,|Tj], taking into account explicitly 
the finite coherence time of the detection process. 
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Figure 2: Real part of the wave packet ( ITH ) at t — 0. The dashed line represents 
envelope ([16]) . 
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